Negative refraction and spatial echo in optical waveguide arrays 
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The special symmetry properties of the discrete nonlinear Schrodinger equation allow a complete revival of the 
initial wavef unction. That is employed in the context of stationary propagation of light in a waveguide array. 
As an inverting system we propose a short array of almost isolated waveguides which cause a relative tt phase 
shift in the neighboring waveguides. By means of numerical simulations of the model equations we demonstrate 
a novel mechanism for the negative refraction of spatial solitons. © 2013 Optical Society of America 
OCIS codes: 190.4370, 190.6135. 
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The optical properties of nano-structured metallic sys- 
tems, often called metamaterials, can be tailored to al- 
most any need [1]. The left-handed materials with neg- 
ative refractive index represent a striking example of 
such kind of artificial structures [2] . It is important that 
the negative refraction is predicted to exist for nonlin- 
ear metamaterials as well [3] . Apparently, the concept of 
negative refraction implies great potential for the control 
and manipulation of optical solitons [1,4]. It should be 
noted, however, that according to the Kramers-Kronig 
relations, the left-handed materials are inherently ab- 
sorptive at the relevant frequencies [5]. Although the 
dielectric nano-structures cannot be described by effec- 
tive medium parameters, it is still possible to apply the 
generalized Snell's refraction law to such systems (e.g. 
see [6]). In this sense, the negative refraction can be re- 
alized in periodic dielectric systems as well [1]. The phys- 
ical mechanisms that govern the electromagnetic fields in 
dielectric nanostructures are different from those in the 
metallic systems, and therefore, the constraint based on 
the Kramers-Kronig relations does not apply. 

Recently, in [7] the negative refraction of optical beams 
at interfaces between bulk dielectrics and tilted photonic 
lattices was studied both theoretically and experimen- 
tally. Moreover, it was shown that the discussed effect 
persists in the presence of nonlinear self-focusing. The 
negative refraction was predicted for linear waves at in- 
terfaces between two photonic crystal waveguide arrays 
as well [8]. Later, in [9] it was demonstrated that the 
solitary wave refraction can be effectively controlled in 
nonlinear optical lattices. Here, we show that the nega- 
tive refraction and spatial echo for solitons can be real- 
ized in nonlinear periodic dielectrics. In particular, the 
soliton propagation in an array of coupled optical waveg- 
uides [4,10] is chosen as a model system. We demonstrate 
that the properly designed scatterer in the system causes 
the desired relative phase shift in the neighboring waveg- 



uides, and leads to the negative refraction of solitons in 
the nonlinear optical lattices. 

Our consideration is based on the fact that in quantum 
mechanics the Schrodinger equation is invariant under 
the simultaneous transformations t ^ —t and t/j ^ tl;"^ . 
Therefore, the conjugation of the system wave function 
at the arbitrary time moment leads to the reversal of the 
system dynamics. Alternatively, it is possible to change 
the sign of the whole Hamiltonian. For example, this 
can be done in the context of cold atoms in optical lat- 
tices by simultaneously changing the signs of hopping 
and interaction constants [11]. The formal analogy of 
the quantum-mechanical systems with light propagation 
in waveguide arrays [12] allows one to expect the similar 
effects in the context of optics. In particular, the com- 
plete image reconstruction in a system of coupled linear 
waveguides was demonstrated theoretically [13] and re- 
alized experimentally [14]. 

In order to explore the above mentioned analogy, let 
us start with the Discrete Nonlinear Schrodinger (DNLS) 
equation: 



dz 



J (^m+l + ^m-l) + 6{z)\^l)m\^ll;r, 



(1) 



where il)m stands for an on-site function of its argument, 
and 0{z) is a stepwise function defined as 6^(2; < 0) = —1 
and 0{z > 0) = 1. It is easy to see that this equa- 
tion is invariant under the simultaneous transformations 
V^m (— Ij^V^m and z ^ —z. Furthermore, let us sup- 
pose that Eq. (1) is valid in the interval |^| > ^ and allow 
the system to evaluate from the initial state iprn{—L) to 
i^ra{—^)' Now, if the system can somehow be brought 
to the state t/jmi^) = (— l)^'^m(— the initial state re- 
emerges at z = L with il)m{L) = {—l)^iljra{—L). The 
transition from il)m{—^) to il)m{^) = {—'^)^i^m{—^) can 
be achieved by evaluating the system according to the 
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following equation: 

dz 



2£ 



(2) 



in the interval \z\ < i.lt should be noted that the neigh- 
boring sites gain tt phase difference in this interval. 

Therefore, for our purpose, the waveguide array de- 
sign must be governed by (1) in the interval \z\ > 
while for \z\ < £ the approximate model equation (2) 
must be valid. The discrete equations (1) and (2) repre- 
sent good approximations for the cold atom dynamics in 
deep harmonic potential or light propagation in optical 
waveguide arrays. A more general theoretical framework 
is given by: 



dz dx^ 



(3) 



for \z\ > £. Here, 6n = Wcos(x), and x(^) = —x{~^) 
stands for the Kerr- type nonlinearity constant. More- 
over, for \z\ < £: 



dz dx^ 



(4) 



where 5n = W cos{x) + (7r/4£) cos(x/2). In these equa- 
tions the variables x and z are scaled in units of 1/i^ 
and kjK^^ respectively. K represents the inverse spac- 
ing of harmonic transverse modulations of the refrac- 
tive index, and k is the carrier wavenumber given by 
k = uquo / c. Here, uo is the laser beam frequency, and 
no is the averaged refractive index. Note that, 5n = 
k'^ [n{x) — no] /noK'^ defines the scaled relative variation 
of the linear refractive index in the corresponding regions 
of the waveguide array. 

Let us first examine Eq. (3). Note that, as far as the 
case of deep lattices is concerned, the wavefunction can 
be assumed to be localized around the maximums of re- 
fractive index profile, i.e. at the points x = 27rm (m is 
integer). Moreover, the fundamental mode for the mth 
potential well can be approximated by: 



- exp {-VW(x - 27rm)Vv^} • (5) 



In this case the tight-binding approximation applies [IS- 
IS], and the full wave function ^(x, z) can be written as 
follows: 

^(X, Z) = ^ llJm{z)(pm{x). (6) 
m 

Furthermore, assuming that the mode overlap of the 
neighboring sites is small, (3) reduces to the DNLS equa- 
tion (1) where the couphng constant J can be easily de- 
termined from the following expression [15]: 



J 



I 



Sn{x)(pm(p- 



Vn+l 



dx dx 



dx. (7) 



Here the normalization of eigenmodes (pm in (5) is chosen 
such that in DNLS equation (see Eq. (1)) the nonlinear 
coefficient value becomes exactly one. 



1 

X=0.002 1 

\- 


1 

7=0 1 1= - 


0.002 


1 1 _ : 


1 
1 
1 


1 
1 
1 


a) 




Fig. 1. Light intensity distributions for a) and b) repre- 
sent the results of numerical simulations with different 
boundary conditions of (3) for L < |z| < ^ and of (4) 
for \z\ < £. Here, L = 750 and £ = 50. The values of 
nonlinear coefficients are shown in the graphs. Dashed 
vertical lines indicate the boundaries of the scatterer. 
c) schematically displays the refractive index profiles in 
the corresponding arrays (see the text for the values of 
modulations). 



In the case of Eq. (4) the lattice of depth is mod- 
ulated with the amplitude tt/M. As far as this modula- 
tion is assumed to be small, similar to Eq. (5) we can 

write (pm{x) ~ exp | — vW^(x — 27rm)^/A/8|- Neverthe- 
less, the harmonic modulation affects the corresponding 
eigenenergies of the modes and leads to a staggered en- 
ergy structure. Using again the expansion (6) for large 
values of W we obtain the negligible overlap between 
neighboring modes J ^ 0. Therefore, we can equate the 
hopping terms to zero which finally gives the discrete 
equation (2). 

We performed the numerical simulations based on 
Eqs. (3) and (4) with the initial conditions at z = —L 
{L £). Initially, the launched spatial soliton propa- 
gates along the trajectory z = vx. Here, the parame- 
ter V determines the direction of the optical beam. Due 
to the effective index mismatch, the scattering process 
takes place at the interfaces of the waveguide array re- 
gion —£< z < £. That causes the change of the beam 
propagation trajectory to z = —vx, and eventually, re- 
sults in the negative refraction of the soliton. Moreover, 
the arbitrary initial profile of the wave packet dit z = —L 
is expected to be reproduced at z = L. 

In the numerical simulations (see Fig. 1) we choose 
L = 750 and ^ = 50. In the region £ < \z\ < L the refrac- 
tive index modulation is Sn = cos(x), and the nonlinear 
coefficient is |x| = 0.002 (for —L < z < —£ we have 
X < 0, while X > foi" ^ < ^ < ^)- Moreover, in the 
range < refractive index is modulated according to 
Sn = 3cos(x) + (tt/M) cos(x/2), and the nonlinear coef- 
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Fig. 2. The same as in Fig. 1 but with two sohtons and 
two scatterers. Note that the nonhnearity is negative in 
the region between the scatterers while it is positive on 
the left and right sides. 

ficient X = 0- The phase modulation of the initial wave 
function determines the spatial soliton propagation "ve- 
locity". As expected, Fig. 1 (a) shows that the soliton 
propagation is not affected by the the thin scatterer if 
the incident angle is normal {v = 0). Furthermore, Fig. 1 
(b) demonstrates that the soliton with v ^ experi- 
ences the negative refraction after the interaction with 
the scatterer. 

Nevertheless, although the output intensity distribu- 
tion is the same among waveguides as at the input, the 
phase distribution is modified. In particular, as it was 
stressed above, t/jmiL) = {—l)^iprn{—L). In order to 
have the exact revival of the wavefunction the second 
scatterer must be introduced. Then, if the second scat- 
terer is identical to the first one, the soliton "lensing" 
can be realized (see Fig. 2). 

It should be stressed that the presented analysis of the 
refraction effects in 1 + 1 dimensions can be extended to 
the case of 1 + 2 dimensions as well. In particular, if 
the propagation direction coincides with the z-axis, the 
scatterer in 1 + 2 dimensions represents a thin film in xy 
plane. The refractive index profile is given by: 

TT 

Sn = W' cos(x) cos(^) + — cos(x/2) cos(?//2), (8) 

and the thickness of the film is 2i. 

In summary, the presented results of numerical simu- 
lations clearly demonstrate that the negative refraction 
of spatial solitons can be realized in an array of cou- 
pled optical waveguides. This effect is related to more 
general phenomenon which can be interpreted as spatial 
echo. Based on the analysis of the theoretical model we 
suggest the inverting devise that allows the complete re- 
vival of the initial wave packet. It must be stressed that 
the discussed effects may have direct implications for the 
telecommunication applications. Similar effects can also 
be realized in different physical systems which can be de- 
scribed by the discrete nonlinear Schrodinger equation. 
These include Bose-Einstein condensates in deep optical 
lattices, coupled defects in photonic crystals, etc. 
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